Simple proof of fault tolerance in the graph-state model 
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We consider the problem of fault tolerance in the graph-state model of quantum computation. 
Using the notion of composable simulations, we provide a simple proof for the existence of an accu- 
racy threshold for graph-state computation by invoking the threshold theorem derived for quantum 
circuit computation. Lower bounds for the threshold in the graph-state model are then obtained 
from known bounds in the circuit model under the same noise process. 



O 
O 

(N 



: 

> 

o . 

m ■ 

cn : 
o . 

o ■ 
^ : 

Ok 
> ' 
G ■ 
d • 

s : 



I. INTRODUCTION 

Alternative models for quantum computation based on 
projective measurements 0, H have recently attracted 
much attention. A common concept of thcse models is 
the simulation of individual quantum circuit operations 
and how simulations can be composed together [j] . More 
specifically, in these models a sequence of single- or two- 
qubit measurements is applied to a collection of fixed ini- 
tial quantum states thereby in effect simulating unitary 
transformations on a smaller subspace of states. 

Although our results are applicable to a much larger 
class of measurement-based models, our analysis will fo- 
cus on a variation of Rausscndorf and Bricgcl's one-way 
quantum computer model pj where computation is per- 
formed by single-qubit projective measurements on some 
initial graph state [5J- Henceforth, we will refer to this 
model as the graph-state model and the computation re- 
alized in it as a graph-state simulation. 

The graph-state model offers a decomposition of a 
quantum algorithm in terms of alternative elementary 
primitives, as well as potential advantages in certain 
physical implementations. For example, suppose entan- 
gling gates can only be realized nondeterministically with 
flagged faults as, e.g., in optical quantum computation 
0. Then, graph-state simulation offers much advantage 
since entangling gates are only used for the preparation 
of the graph state, which can be done independently from 
the main computation 0. 

In any physical realization of quantum computation, 
unknown errors will always be present and they will have 
to be corrected using quantum error-correcting codes in 
a fault-tolerant manner. An important question is there- 
fore under what conditions computation can be executed 
reliably in the graph-state model in the presence of phys- 
ical noise. In the circuit model, fault-tolerant methods 
; 8j] are available for the reliable execution of any desired 
computation, if the noise is sufficiently weak. Now, if 
such a fault-tolerant circuit is simulated in the graph- 
state model with sufficiently weak noise, will the same 
desired computation be reliably executed? More specifi- 
cally, this is answered by first analyzing noisy simulations 
of individual operations and then how the noisy simula- 
tions compose together. 

The first results on this problem were reported in the 
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Ph.D. thesis of Raussendorf IQ]. This work proved the 
existence of an accuracy threshold for cluster-state com- 
putation for various independent stochastic error mod- 
els. More recently, Nielsen and Dawson E3 obtained 
proofs for the existence of an accuracy threshold in the 
graph-state model that apply to more general error mod- 
els (including errors due to nondeterministic gates) by re- 
duction to a threshold theorem for local non-Markovian 
noise 11}. In addition, they established a conceptual 
framework and two technical theorems that are of inde- 
pendent interest. 

In this paper, we use the concept of composable simula- 
tions 
model 

the existence of an accuracy threshold in the graph-state 
model. Furthermore, for any specific form of noise, our 
proof allows known lower bounds on the threshold in the 
circuit model to be translated to equivalent bounds in 
the graph-state model. We discuss in particular how the 
two lower bounds are related for commonly used fault- 
tolerant architectures based on self-dual CSS codes. 



II. REVIEW OF THE GRAPH-STATE MODEL 

We begin by briefly introducing the graph-state model 
in the ideal noiseless case. Various recent interpretations 
of this model have been reported and reviewed recently 
0,JP, 0, [H, EÍ, Hj- We follow the language in 
Ref. [4J , which explicitly uses the notion of composable 
simulations that forms the core of our subsequent analy- 
sis. Our discussion in this section is intended to introducc 
the bàsic notions and terminology that we will use later 
in our proof. 

In the circuit model, an arbitrary quantum computa- 
tion can be decomposed into state preparation, measure- 
ments, and a universal set of gates. To show the univer- 
sality of the graph-state model of quantum computation, 
it suffices to show that (i) each element for universality in 
the Standard model can be simulated and (ii) the simula- 
tion can be composed to simulate the entire computation. 
The approach is to first define an appropriate notion of 
simulation that is composable, followed by a complete 
recipe to composably simulate each element needed for 
universality. 



2 



We first describe the notion of composable simulations. 
Lct T be an operation (a superoperator, or completely 
positive trace-preserving map) to be simulated and 
be the associated operation that simulates T . For sim- 
plicity, let T act on n qubits. In the general case, T can 
have quantum and classical input and output of arbitrary 
dimensions, but this only requires extra notations and 
therefore will not be written out explicitly here. For a 2n- 
bit string x, let V x be the superoperator corresponding 
to conjugation by the Pauli operator indexed by x. Our 
composable simulation Sjr takes two inputs, a classical 
2n-bit string e- m and an n-qubit quantum state V eiIl (pi n ), 
so that V/9i n , Vein, it acts as 

PzaMt e OUt ®(^e out O^)(p in ), (1) 



where e ou t is some new 2n-bit string that appears with 
probability p £mt ■ (Throughout the paper, the symbol for 
a bit string such as e- ln also labels the corresponding den- 
sity matrix.) To rephrase the above definition, for each 
specific classical output e out , <Sy evolves the arbitrary 
state pi n according to the intended operation T up to a 
new known succecding Pauli operation V eout , despite the 
Ve ln oceurring to p- m prior to the simulation. 

Note that e ou t is a function of e; n and the measurement 
outeomes obtained in Sf, and this function depends on 
iSjr. However, the statistics of e ou t has no consequence, 
because composable simulations work for all measure- 
ment outeomes and for all ei n — all outeome histories lead 
to vàlid simulations, where an "outeome history" denotes 
the set of all measurement outeomes collected in a spe- 
cific run of the simulation. As we will see next, this is 
important as it will allow us to compose simulations of 
individual operations to obtain a simulation of the com- 
bined operation. 

Now, consider simulating a sequence of l operations 
{J-j}, and we will see that it can be done by composing 
the sequence of simulations {<Sj^}. By repeated applica- 
tions of Eq. QJ, Vpi n , Vei n , 



5j7, o • ■ • o S ri (e in íg) V eitl (Pm)) 

eout 



(2) 



which states that, for all outeome histories, the entire 
sequence of operations {J-j} is simulated properly, up to 
a final overall V eout (which just redefines the final classical 
outeome of the computation). 

We will now describe how composable simulations are 
realized in the graph-state model. Let T denote a graph 
with vèrtex set V^r) and edge set E(T). One way to 
specify and to create the graph state corresponding to 
T is to start with the initial state ®iev(r) 1+) aim then 
apply a controlled-phase (cphase) gate to each pair of 
qubits in E(T) (where CPHASE [aò) = (-l) ah |aò) in the 
computation basis). In other words, each vèrtex corre- 
sponds to a qubit initially in the state |+), and each edge 
corresponds to a subsequent CPHASE. 
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As precursor to a graph-state simulation, our next step 
is to composably simulate a universal set of circuit ele- 
ments (state preparation, measurements, and a univer- 
sal set of gates), using single-qubit measurements and 
CPHASE. In the circuit model, it suffices to prepare any 
Pauli eigenvector and measure along any Pauli basis. 
Both of these can be trivially simulated in the graph- 
state model using single-qubit measurements. For the 
universal set of gates, we choice the Clifford group gen 
erators {H, S = e 
non-Clifford T = e 
Standard Pauli operators. Figurenshows how to com- 
posably simulate these gates, with the classical registers 
omitted for simplicity. In Fig. ^ qubits are represented 
as circles. The boxed circles contain the quantum inputs, 
unboxed ones are prepared in |+), and open circles (un- 
measured qubits) contain the quantum outputs. Edges 
denote cphase gates acting on the adjoined qubits. The 
measurement bases for each qubit are given in the cir- 
cle. The quantum state at the input of each pattern has 
known Pauli corrections labeled by the classical regis- 
ter ei n (not shown), which depends on past measurement 
outeomes. In the simulation of T, e; n is used to control 
one of the quantum measurements. The output quantum 
state also has Pauli corrections labeled by an updated 
string e ut- Each simulation pattern defines an update 
rule, mapping ei n and measurement outeomes obtained 
in the pattern to e ou t. 
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FIG. 1: Composable simulations for (a) the Hadamard gate (H), 
(b) the rotation around the 2-axis by tt/2 (S), (c) the CPHASE and 
(d) the rotation around the 2-axis by 7r/4 (T). Note that we use 
the CPHASE to simulate itself, since it can be built in as a vertical 
edge of the graph. We have omitted the input classical registers 
and their updates for simplicity. The symbols Mx and My indicate 
measurements of <r x or <r y on the corresponding qubits, and Mt 
indicates a measurement of the observable (o- x ±<Ty)/\/2 depending 
on whether there is a <r x correction in the input qubit. 



Any circuit (sequence of gates and measurements on 
Standard initial states) can then be simulated by compos- 
ing a sequence of simulations by identifying the quantum 
output of one simulation (the open circle) with the in- 
put to the next (the boxed circles) and similarly for the 
classical registers. The combined simulation thus con- 
sists of single-qubit measurements on qubits prepared 
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in a graph state (with the CPHASE being part of the 
graph state preparation) , giving a complete recipe for 
the entire graph-state simulation. Note that evolutions 
of single qubits and their interactions (cphase) in the 
simulated circuit are represented in the graph as linear 
paths and the links between them, respectively. As an 
example, Fig.[21 shows how a composition of the mea- 
surement patterns for the simulation of H and cphase 
leads to a new pattern for the simulation of the operation 

CNOT = (I <g) íf)CPHASE(/ ® H). 



Ah 



FIG. 2: A schematic diagram of the composition of the pat- 
terns in Fig.^a), Fig.^c), and Fig.^a) that simulates CNOT = 
(7 (g> H) CPHASE (/ (g> H). In the dashed ellipse on the left, the out- 
put of the measurement pattern simulating the first H is identified 
with the input for the lower qubit of the CPHASE simulation, whose 
output for the same qubit is identified with the input qubit of the 
simulation of the second H (right ellipse). On the right is the result 
of the composition. 



III. NOISY GRAPH-STATE COMPUTATION 

We now investigate how noise at the level of the graph- 
state simulation maps to noise in the simulated oper- 
ations and, most importantly, whether such "simulated 
noise" can be tolerated by simulating a fault-tolerant cir- 
cuit described in the circuit model. We begin by men- 
tioning a modification to the graph-state model that is 
necessary for fault tolerance. Since the ability to prepare 
fresh qubits and interact them with the existing ones is 
essential for all fault-tolerant constructions |2Ï|, instead 
of creating the entire graph state before computation, 
a minimal modification to the model is to build the re- 
quired graph state dynamically as the computation pro- 
ceeds 0, 0, Oa, 113 • The simulated circuit defines a partial 
time ordering of the simulations and the measurements 
used therein, inducing a partial ordering of the qubits in 
the graph state. The qubits can be added slightly be- 
fore their preceding neighbors are measured, as long as 
the CPHASE gates are applied according to the time or- 
dering of the simulations. This change in the model still 
preserves the appealing feature of the graph-state model 
in that all unitary interactions are applied prior to and 
independent of the measurements that realize the com- 
putation. 

Corning to the main part of this paper, we must an- 
alyze how physical noise affects the elementwise simula- 
tions and how the noisy simulations compose together. 
The elementary steps in the simulation are the prepara- 
tion of |+), the cphase, the single-qubit measurements, 



and the storage of qubits. Moreover, each operation be- 
longs to a unique simulation. Thus, noise afiheting a 
given operation only acts within one simulation. In par- 
ticular, an erroneous cphase cannot affect two successive 
simulations. 

In any noise model and without loss of generality, each 
noisy state preparation or noisy gate can be expressed as 
the ideal operation followed by a noise operation. Hence, 
noise operations intersperse pairs of successive ideal oper- 
ations. A noise operation is a system-environment cou- 
pling, and it can always be described by some unitary 
joint evolution 



t^fault = 



(3) 



where Pi ranges over all nontrivial Pauli operators in- 
dexed by i acting on the output system of the preceding 
ideal operation and each Ai is an arbitrary operator act- 
ing on the environment, subject to the condition that 
C^fauit is unitary. A noisy measurement is modeled as the 
ideal measurement preceded by a noise operation given 
by Eq. Q, with Pi acting on the qubits to be measured. 

We first consider independent stochastic noise pro- 
cesses. In this case, each noise operation is by assump- 
tion acting on a separate environmental register, which is 
mapped to orthogonal states by the two terms in Eq. © . 
Physically, this assumption corresponds to the require- 
ment that a record be kept in the environment whenever 
faults oceur, which can in principle be read to indicate 
the location of faults. In more detail, the two terms re- 
sult in perfectly distinguishable environmental states, so 
that the corresponding states in the system do not in- 
terfere with one another, and their normalization can be 
interpreted as the probabilities of the first or second term 
in Eq. oceurring. These two terms thus correspond to 
the two events of not having or having a fault. We call 
the second term in Eq. the fault operator or simply 
the fault. A fault path for the entire computation is an 
event oceurring with some definite probability describing 
whether each noise operation results in a fault or not. 

Our first goal is to show that faults within one simu- 
lation only affect that simulated operation, even though 
classical registers that carry the Pauli corrections and 
control the simulation are shared by many simulations. 
Consider a sequence of simulations {Sy^ j } applied to an 
input J2 e Pe in Gín <8> Pe ia (pin)- Suppose some number of 
faults oceur within . Each term in the expansion of 
Eq. © of all these fault operators consists of Pauli oper- 
ators acting on the simulation qubits which can be com- 
muted to the end of the simulation (since, as shown in 
Fig.n each simulation is realized by a sequence of uni- 
tary CPHASe(s) and single-qubit measurements). This 
results in a combined fault operator, each term in the 
Pauli expansion of which contains some Pauli operator 
acting on either the output classical registers of or its 
quantum output, or both. The most general erroneous 
output is thus given by J2e mt pioL e out ® Amt for some 
distribution {p22 t }, where e ou t is some possibly erro- 
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neous classical output and p ou t = ^t^e 1 " °·^ r i(Pm)), 
£ Gout is the completely positive trace non-increasing map 
induced by the combined fault operator on the quan- 
tum output and is conditioned on e ou t, and eò„ e t al la- 
bels the ideal corrections at the output in the absence 
of faults inside Sp 1 and depends on ej n . Let £ eout = 
^eout ° ^e out ° V e ideai. Then, the output of can 

be rewritten as J2 ecmi , Pelt e ou t ® ^,(4». ° -^(Pm))- 
Hence, besides the extra £ eout , the noisy output state is 
of the same form as some ideal noiseless output, with the 
classical register reflecting the Pauli correction on the 
quantum state. In particular, this means that we can in- 
clude errors in both the quantum and classical registers 
in £ Eout and intèrpret it as a simulated fault operation 
following the simulated T\. The above analysis can now 
be repeated to subsequent simulations, so that a simu- 
lated fault appears after each erroneous simulation. In 
each term labeled by e ou t , the simulated evolution on the 
system and the environment is the intended computation 
(the sequence {J~j}) interspersed by the action of simu- 
lated fault operators (whose particular type may depend 
on ei n or e ou t at the corresponding erroneous simulation) . 

We pause to discuss the above argument again. The 
composable simulation has been described in many dif- 
ferent ways in the literature, such as feed-forward of mea- 
surement outcomes and propagation of by-product Pauli 
operations. Since the classical knowledge (correct or not) 
of these by-product Pauli operations from one simulation 
step is input to the next, it is worrisome that an error in 
them will feed forward, inducing highly correlated simu- 
lated faults in the simulated circuit, even if initial faults 
in the simulation are uncorrelated. It only takes a shift 
in one's perspective and inspection of the composability 
rcquirement to recognize a simpler interpretation of the 
error action. In particular, errors in the classical informa- 
tion of the by-product Pauli operations e ou t are equiva- 
lent to unknown Pauli errors in the quantum output p ou t 
of the erroneous simulation. The above argument takes 
full advantage of the cquivalcncc and mathematically re- 
defines e ou t to indicate the by-product Pauli operation, 
attributing any "mismatch" with an ideal simulation to 
noise acting on the quantum output of the simulation of 
T\ alone. From this point of view, the errors in classical 
information are localized and do not propagate. Being 
able to localize errors to individual simulated operations 
achieves a simple and direct mapping from the noise in 
the simulation to noise in the simulated circuit. 

We can now finish the proof of the existence of an 
accuracy threshold for independent stochastic noise in 
the graph-state model, using the threshold theorem for 
Standard quantum computation 0, 0, 0, 0] : In the 
circuit-model proof, certain fault paths are "good" and 
can be proved to give the ideal computation results. 
"Bad" fault paths form the complement of the "good" 
ones and have suppressed probability if the physical fault 
probability is below a certain critical value, the accuracy 
threshold. Consider the noisy graph-state simulation of 



a fault-tolerant circuit. In the final output of the fault- 
tolerant circuit simulation, consider each e ou t term. Our 
arguments based on composability ensures that the evo- 
lution of the quantum state is simply the intended sim- 
ulated operations, interspersed by the action of faults. 
Since each fault path in the simulation is mapped to a 
unique fault path in the simulated circuit due to error 
localization, good fault paths in a graph-state simulation 
can be defined as those resulting in good fault paths in 
the simulated circuit j^. All other fault paths in the 
simulation are bad, and their probability will be sup- 
pressed below a certain accuracy threshold just as in the 
circuit model, because a simulated fault appears after 
some simulated operation only if there is at least one 
fault in its simulation. Furthermore, the probability of 
this happening is at most the sum of the fault probabili- 
ties of all the elementary steps in the simulation. Then, 
with reference to Fig. ^ we note that the simulation of 
each gate in our universal set involves the use of one to 
two CPHASE gates and zero to two measurements. Thcre- 
fore the probability of any simulation containing faults is 
bounded by p s im < 4p. More specifically, if po is the 
threshold value of the fault-tolerant architecture used in 
the circuit model and if p < po /4 in the simulation, then 
Psim < Po and the final measurement outeome will pro- 
vide the correct computation results with the desired ac- 
curacy. This holds for each e ou t term in the final state 
of the simulation, thereby establishing a threshold lower 
bound of po /4 for the graph-state model. 

In the above, we have related the accuracy threshold in 
the graph-state model to that in the circuit model by the 
direct simulation of fault-tolerant architectures designed 
in the latter. However, we note that, in order to obtain 
the above threshold bound, we assumed that the fault- 
tolerant simulated circuit makes use of the same universal 
set as ours. In general, the same gate sets need not be 
used in both models, and elementary measurement pat- 
terns need to be composed to simulate a single operation 
in the simulated circuit. In particular, in most studies, 
CNOT rather than CPHASE is used as the elementary in- 
teraction. In this case, the measurement pattern in Fig.[21 
for the simulation of CNOT implies the threshold condi- 
tion p < po/5. However, in many cases of interest this 
lower bound is pessimistic. For example, in fault-tolerant 
designs based on self-dual CSS codes (e.g., [Ï5ll23ll2^|). 
CPHASE can replace CNOT as an alternative bitwise en- 
coded operation and can also be used in error correction 
with a small number of additional H gates. Since there 
is no overhead for simulating single-qubit state prepa- 
ration, measurement, or the CPHASE in the graph-state 
model, the thresholds for circuits based on these codes in 
the circuit and graph-state models will be essentially the 
same. 

We now proceed to prové the existence of an accuracy 
threshold for the graph-state simulation for local non- 
Markovian noise. We will make use of our observation of 
the localization of errors and the threshold results in the 
circuit model [TTL Íl5| ■ 
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In the local non-Markovian error model, the noise op- 
erations still have the form given by Eq. © and they act 
on the system in the same way as in the local Markovian 
model. However, different noise operations may now act 
on the same environmental register, and the term acting 
trivially and nontrivially on the system may not map the 
environmental register to orthogonal states. Altogether, 
faults can combine coherently. Furthermore, a fault no 
longer corresponds to an "event," in the sense that proba- 
bilities cannot be assigned. Instead, one imposes that the 
strength of the fault operator at each location is boundcd 
below a certain value r\ — i.e. 1 1 £\ Pi ® Ai\ Isup < V- 

To simplify the analysis, we consider the purification 
of the graph-state simulation, where measurements are 
replaced by coherent operations by attaching extra an- 
cillary qubits. In our noise model, noisy measurements 
are modeled as being ideal with noise factored into the 
preceding noise operations, so that changing our descrip- 
tion of the measurements does not affect the analysis. 
Likewise, the classical 2n-bit string carrying e- m can be 
replaced by a 2n-qubit register in the state \e- m ) and any 
adaptive operations inside these equivalent simulations 
will be controlled by this quantum register. The update 
of this register to obtain |e ou t) can also be done coher- 
ently by controlling gates from the extra ancillary qubits 
and also by doing the classical processing reversibly. We 
emphasize that this alternative coherent description is 
purely mathematical and is also employed in the circuit- 
model proofs in Refs. [ïïl lïfij . The composable simu- 
lation 5jf is now a conjugation by a unitary operator 
Sp taking two inputs |ei n ) and P ein \ip- m ) and some ancil- 
lary qubits starting in the fixed state |+)® fe j2í|, so that 
VIV'in)) V|ej n ) it acts as 

^(|e ln )®P ei J^ n )® |+)» fe ) 

(4) 

= Ec. |e ut> ® |i) ® |&) <8 P eout F|Vi„) , 

i 

where {|0i}} is the orthonormal basis on which measure- 
ments are to be performed, {\i}} is the computation basis 
with i labeling the possible measurement outcomes car- 
ried by the extra ancillas we have introduced, Cj is the 
amplitude of the ith term, |e ou t) is a 2n-qubit state that 
depends on e- ln and i, and F is the simulated unitary 
operator. 

Having expressed the fault-tolerant circuit to be simu- 
lated as well as the graph-state simulation itself unitarily, 
a unitary noise operation of the form of Eq. J2J is inserted 
at every location in the simulation (where locations are 
specified by the original graph-state simulation before 
the unitary idealization) . The output state is a linear 
superposition of states, each evolved according to a spe- 
cific set of fault operators and expanded in the eigenbasis 
of all measured operators (including both measurements 
part of the graph-state simulation and also measurements 
originally in the simulated circuit). Fault paths can again 
be "good" or "bad," defined as in our discussion for in- 
dependent stochastic noise. For each term evolved by a 
good fault path, a final quantum state that will provide 



the correct statistics will be generated, independent of 
the state of the register \i) coherently carrying the mea- 
surement information due to the localization of errors. 
This is because, for each term in the Pauli expansion of 
faults acting on |e ou t), the register |e ou t) can always be 
taken to carry the correct Pauli correction by redefining 
the error acting on P eout _F|-0; n ) exactly as in our previous 
discussion. Therefore, for each term in this Pauli expan- 
sion, good fault paths in the simulation are mapped to 
good fault paths in the simulated circuit that produce 
the ideal computation results, using the threshold theo- 
rem in the circuit model. Hence, by linearity, the wholc 
coherent sum of these terms will also produce the ideal 
computation results. 

It remains to bound the sup norm of the bad fault 
paths of the graph-state simulation, which can combine 
coherently. Following the threshold theorem in the cir- 
cuit model for local non-Markovian noise 0, ^| , it suf- 
fices to bound the sup norm of the "bad" part of a given 
simulation (i.e., the sum over terms of the form ^ Pi® Ai 
in at least one location within this simulation) . But this 
sup norm is simply bounded by »7sim < 4?y, where rj is a 
bound on the sup norm of the fault operator acting on 
each location in the simulation (by the triangular inequal- 
ity of the sum norm). Thus r\ < 770 /4 is the threshold con- 
dition for the graph-state model if 770 is the established 
threshold strength for the circuit model. 

IV. CONCLUSION 

To conclude, we have invoked the composability prop- 
erty of simulations in the graph-state model to show that 
faults in the graph-state simulation of any quantum cir- 
cuit (and of a fault-tolerant circuit, in particular) can be 
viewed as affecting the simulated operations alone. Thus, 
the existence of an accuracy threshold for the graph- 
state model follows from the threshold theorem in the 
circuit model for the same noise process. As an aside, the 
same insight can be applied to other measurement-based 
models of quantum computation and the teleportation of 
gates. Although other proofs for the existence of an ac- 
curacy threshold in the graph-state model have already 
been reported for a variety of error models 0, , we be- 
lieve our analysis provides an alternative, conceptually 
different and in many respects simpler way of thinking 
about fault-tolerant circuit simulations. 

We note that in optical implementations of graph-state 
computation Q , gate nondeterminism and photon losses 
give additional sources of faults not treated in this work. 
The works in Refs. 261 show how to control these 

faults by preparing microclusters. A precise threshold 
analysis in this setting is pursued elsewhere |27l |. 
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